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The coherence length of superflnid neutron matter is calculated from the microscopic BCS wave- 
function of a Cooper pair in momentum space making use of the Bonn meson-exchange potential. 
We find that the coherence length is proportional to the Fermi momentum-to pairing gap ratio, 
in good agreement with simple estimates used in the literature, and we establish the appropri¬ 
ate fitting constants using our numerical data. Our calculations can be applied to the problem of 
inhomogeneous superfluidity of hadronic matter in the crust of a neutron star. 
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Calculations based on BCS theory with phenomeno¬ 
logical nucleon-nucleon (NN) forces indicate that neu¬ 
tron matter is superfluid in a wide region of densi¬ 
ties and temperatures. In particular, at the densities 
corresponding to the interior of a neutron star crust 
(4 X 10^^ gcm“^ < p < 10^"^ gcm“^) neutrons couple 
in the singlet isotropic channel while at larger den¬ 
sities they pair in the ^P 2 state 1^. Although many in¬ 
vestigations have been devoted to the superfluidity and 
superconductivity of neutron, /3-stable and nuclear mat¬ 
ter, little attention has been paid to possible effects of 
inhomogeneities in hadronic superfluids or superconduc¬ 
tors. In fact the interior of neutron stars, the subject of 
most of such calculations, is often referred to as the only 
existing example in nature of infinite superfluid neutron 
or /3-stable matter, since in an atomic nucleus the wave- 
functions of the Cooper pairs are limited in extension by 
the potential well. On the other hand, superfluidity in 
a neutron star crust represents a case intermediate be¬ 
tween the nucleus and the idealized infinite system, since 
superfluid neutrons in the inner crust occupy a region 
where a lattice of nuclei creates strong inhomogeneities 
in the medium. From the point of view of astrophysi- 
cal observations, probably most of the visible effects of 
the presence of superfluids in a neutron star are due to 
phenomena in the crust ^ . 

An important length scale of the neutron superfluid is 
the coherence length. From a microscopic point of view 
the coherence length represents the squared mean dis¬ 
tance of two paired particles (a Cooper pair of neutrons) 
on top of the Fermi surface. The magnitude of this quan¬ 
tity affects several of the physical properties of a neutron 
star crust. First of all, neutrons paired in a singlet state 
form quantized vortices induced by the rotational state of 
the star. These can pin to the nuclei present in the crust, 
possibly leading to the observed sudden release of angu¬ 
lar momentum known as pulsar glitches. The magnitude 
of the pinning force depends on the size of the vortex 
cores, which is equal to the coherence length of the neu¬ 
tron superfluid. A second question is how properties of 


the neutron superfluid change due to the inhomogeneous 
environment of a neutron star crust, a problem related to 
the average thermodynamical property of neutron mat¬ 
ter 1^,^. In the inner crust, depending on density, nuclei 
of different shapes and sizes are present. At a density of 
~ 10^^ gcm“^, spherical nuclei cease to be energetically 
favored and are replaced first by cylindrical nuclei, then 
slabs to end up with holes, where the roles of protons 
and neutrons are exchanged, see Refs. [||Q| for further 
details. Only at higher densities, corresponding to what 
is called the core of the star, do nuclei merge into the 
uniform medium. The fact that neutron superfluidity in 
neutron star matter is actually a problem of inhomoge¬ 
neous superfluidity in hadronic matter has been noticed 
quite recently ||,p. According to Anderson’s theorem 
[H the electron density of states in a superconductor is 
changed very little from a pure metal to an alloy of sim¬ 
ilar chemical properties. The physical situation we are 
examining here is quite different, since both the average 
density of states and the effective neutron-neutron ma¬ 
trix elements are changed compared to the uniform case 
when one considers the presence of nuclei. 

The typical dimension of nuclei in the inner crust of a 
neutron star is Rn ss 4 — 6 fm. This number is, in an ap¬ 
propriate range of densities, comparable to the coherence 
length ^ as estimated from existing BCS calculations. If 
Anderson’s theorem holds also for neutron star matter, 
(limit where Rn <C there will be no appreciable vari¬ 
ations in the superfluid properties induced by the nuclei. 
On the other hand, if Rn ^ superfluid will change 
its properties locally. This limit has been investigated in 
a recent series of papers Q where it was found that some 
thermodynamical properties like e.g. the neutron specific 
heat may change by a very large amount. Unfortunately, 
the situation is complicated by the fact that Rn and ^ 
are of the same order of magnitude. 

Clearly, the coherence length represents a critical pa¬ 
rameter by which one can establish the behavior of an in¬ 
homogeneous superfluid. It sets the scale for the possible 
spatial variation of the pairing properties of the system. 
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and thus plays a role if some inhomogeneities are present 
in the system at a length scale comparable to it. 

A simple estimate of the coherence length is obtained 
by assuming constant matrix elements between particle 
states within a shell centered at the Fermi momentum 
and zero outside (see for example Ref. [||). This gives 





( 1 ) 


where kp is the Fermi momentum and is the pairing 
gap at the Fermi momentum. The value of the constant 
K depends on the approximations used, and different 
values have been reported in the literature. However, 
the estimate (|^ does not consider thoroughly the matrix 
elements of the particle-particle interactions, and in prin¬ 
ciple can be too rough if one wants to calculate the co¬ 
herence length at all values of the density. In view of the 
important astrophysical considerations discussed above, 
it is desirable to calculate the coherence length from a 
microscopic study of the wavefunction of a Cooper pair. 
The aim of this Brief Report is thus twofold: firstly we 
check the validity of the simple estimate and search 
for possible deviations or fluctuations from its average 
behavior. Secondly, in those cases where (|^) is well re¬ 
produced, we determine the best value of K and check 
the dependence of that formula on the neutron interac¬ 
tion in the particle-particle channel. To this end, we shall 
make use of the Bonn A potential. This meson-exchange 
potential is defined by the parameters of Table A.l in 
Ref. @. 

Our microscopic calculations of coherence lengths have 
been carried out as follows. Let (j){v) be the wavefunction 
of the relative motion of the two neutrons in a Cooper 
pair, r being the relative coordinate of the two particles. 
The coherence length ^ is given by 


2 ^ J d^x\^{r)\^r^ ^ 

f (Px\(j){Y)\'^ f d^k\x{^)\'^ 
/ q °° dkk'^\dxik)/dk\'^ 


where in the second equation the expectation value is 
calculated in momentum space, x(k) being the Fourier 
transform of Eq. (||) is particularly suited for nu¬ 

merical computation, since the BCS equations for a uni¬ 
form system are solved in momentum space. The wave- 
function of the Cooper pair in momentum space is given 
by (apart from an unimportant normalization constant) 


X{k) = 


A(fc) 

E{k) 


( 3 ) 


where A(A;) is the fc-dependent pairing gap, while E{k) = 
A/A2(fc)-p(e(fc)-^)2 is the quasi-particle energy and 
e(fc), jj, are the neutron single-particle energy and chemi¬ 
cal potential respectively. Notice that due to the isotropic 
nature of ^S'o pairing all quantities depend only on the 
magnitude of the momentum. 


We calculate the pairing gaps and wavefunctions in 
momentum space as functions of the momentum from 
the BCS gap equation [|] 


A{p) = 


1 

2 


d^ k 

(^ 


v{p, k) 


A(fe) 

E{k) 


( 4 ) 


for different values of the Fermi momentum. Here v{p, k) 
are the matrix elements of the NN interaction, ob¬ 
tained from the Bonn A potential, following the proce¬ 
dure outlined in Ref. |^. The single-particle energies 
e{k) were obtained from Brueckner-Hartree-Fock calcu¬ 
lations with the same potential, see Ref. for details. 
The chemical potential p, was set equal to the Fermi en¬ 
ergy e{kF). Notice that for our purposes the pairing gap 
has to be calculated at all wavenumbers and not only 
at the Fermi surface, as in the simple formula (P. The 
wavefunctions x(^) momentum space for four differ¬ 
ent values of the Fermi momentum {kf =0.1, 0.4, 0.8 
and 1.2 fm“^) kp are shown in Fig. 1. The square of the 
derivatives of the respective wavefunctions are displayed 
in Fig. 2. It is seen that the spread of the wavefunc¬ 
tion in momentum space < > varies considerably ac¬ 

cording to the value of the Fermi wavenumber kp. The 
distribution in momentum space is particularly wide for 
intermediate values of kp and, since we have through the 
uncertainty relation that <k^>=^^«l/ <r^>,a 
smaller value of the coherence length for these values is 
expected. Due to the limited, but significant variation 
of the the pairing gap A{k) within a single peak (each 
peak occurs at the Fermi momentum), the wavefunction 
is distorted with respect to a simple gaussian and there¬ 
fore deviations from the simple relation (|^ for a gaussian 
distribution should be expected. This becomes evident 
by looking at the squared derivative of the wavefunction 
{dx{k)/dk)^ which shows two unequal peaks. 

In Fig. 3 we present the squared wavefunctions 
r2|(^(r)|2 for the same values of the Fermi momentum 
as in Figs. 1 and 2. The wavefunctions oscillate with a 
wavelength of the order kp^ and extend for quite long dis¬ 
tances compared to the range of nuclear forces (a few fer- 
mis at most). Not surprisingly, this behavior is due to the 
relatively small value of the energy of the pair compared 
to the Fermi energy (a comparison can be made with 
another weakly bound state with a very spread wave- 
function, the deuteron |^). 

In Fig. 4 we show our microscopic calculations of the 
coherence length for pure neutron matter performed us¬ 
ing Eqs. ti- The microscopically calculated points, 
given by the solid line, can be fitted by a function linear 
in the parameter kp/Ap, where Ap is the pairing gap at 
the Fermi surface. We find as a best fit, see Ref. for 
further details on the form for 


, kp 

? = 1 — X 

2m Ap 


( 5 ) 
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where x = m/m* is the ratio between the bare and the ef¬ 
fective mass, and the dimensionless parameter 9 is found 
by a direct fit to our data to be d = 0.814 (for a; = 1) and 
9 — 0.836 (with the correct density dependent value of x 
from the miscroscopic calculation). The choice a; = 1 is 
commonly used in the literature. 

The simple linear behavior of the momentum-to gap 
ratio is in very good agreement with Eq. (|^) and shows 
that a Cooper pair represents a packet having a nearly 
well-defined relative momentum. As seen from Fig. 4, 
the use of a; = 1 (effective mass equal to bare mass) does 
not change the picture very much, as is evident from 
the fact that at e.g., the relatively high Fermi momen¬ 
tum of kp = 1.0 fm“^ the effective mass differs from 
the bare one by less than 4%. However, deviations from 
Eq. d^) are found, as can be seen in Fig. 4, especially 
at higher densities. In Fig. 4 we have also plotted the 
values of the coherence length from an approximation of¬ 
ten used in neutron star studies, namely the choice of a 
coefficient 9 = 2/^ = 0.637 in Eq. (||) [||. Our values 
were on the average more than 20% larger than the re¬ 
sults obtained with this choice. It is therefore interesting 
to see that our microscopic calculation of the gap and 
single-particle properties through a complicated many- 
body scheme, yields a qualitative similar result as that 
of Eq. (1), derived originally in a solid state context |^. 

A relevant point is whether the value of the coeffi¬ 
cient 9 depends on the interaction chosen. Besides, it 
might be that the good agreement with Eq. (^ gets worse 
with other interactions. To partially answer these ques¬ 
tions, we repeated the same calculations making use of 
an effective force, the Gogny force, which pairing prop¬ 
erties have been investigated in several works . The 
Gogny force is an effective interaction fitted to repro¬ 
duce various nuclear data, and contains therefore effects 
of nucleon-nucleon correlations not included in a bare 
nucleon-nucleon interaction. It is therefore both qualita¬ 
tively and quantitatively different from a bare NN inter¬ 
action like the Bonn potential. We find that with the D1 
parameterization of the Gogny force a very good fit can 
be made for Eq. (3) with a coefficient 9 = 0.815, very 
close to the value found with Bonn A. This near-equality 
of the coefficient 9 for the two different interactions is not 
predictable on the basis of simple arguments. However, it 
ought to be stressed that we have not included so-called 
polarization effects in the calculations of the pairing gaps 
of Ref. dp. Such effects are expected to reduce by at 
least a factor of two the value of the pairing gap, see Ref. 
[pd| . How these many-body effects will change the coher¬ 
ence length remains to be investigated. However, using 


the weak-coupling approach of Eq. (|^) for the coherence 
length, one sees that ^ may be twice as large. 

To conclude, we calculated microscopically the coher¬ 
ence length of superfluid neutron matter. Eq. (i with 
the nearly interaction-independent coefficient 9 — 0.83 
represents the main result of our calculations. The data 
reported in this brief report may be useful in quantita¬ 
tive calculations of the superfluid properties of neutron 
star crusts. Moreover, as can be seen from Fig. 4, the co¬ 
herence length is for several Fermi momenta larger than 
the typical size of known nuclei. This poses several con¬ 
straints on the use of local density approximations in the 
study of neutron star crust properties [^]. 
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FIG. 1. Wavefunctions in momentum space, x{k) for four different values of the Fermi momentum, kp ~ 0.1 fm”^, kp — 0.4 
i“^, kp — 0.8 fm“^ and kp = 1.2 fm~^. The wavefunctions peak at the corresponding value of the Fermi momentum. 
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FIG. 2. Plot of \dx{k)/dk'\^ proportional to the spread in space of the Gooper pair. Legend as in Fig. 1. 
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FIG. 4. Coherence length for neutron matter calculated numerically with the Bonn potential (solid line) using Eqs. (2)-(4). 
The line with long dashes are for a ^ which fits the numerical points with 0 — 0.814 while for the short-dashed line 0 — 0.836, 
see text for further details. The dotted line represents ^ given by an approximation frequently used in the vortex pinning 
literature [9]. 
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